


C/3 ■ 

o ■ 

• I— I ■ 
-(— > . 

CI.. 

9: 

^ '. 

o ■ 

•I— I , 
^ , 

Oh- 



> 

m 

o 



X 



Quasi-phase-matching of high harmonic generation using polarization beating in 

optical waveguides 

Lewis Z. Liuo Kevin O'Keeffe, and Simon M. Hooker 

Clarendon Laboratory, University of Oxford Physics Department, 

Parks Road, Oxford 0X1 3PU, United Kingdom 

(Dated: April 17, 2013) 

A new scheme for quasi-phase matching high-harmonic generation is proposed in which polariza- 
tion beating within a hollow core birefringent waveguide modulates the generation of harmonics. The 
evolution of the polarization of a laser pulse propagating in a birefringent waveguide is calculated 
and is shown to periodically modulate the harmonic generation process. The optimum conditions 
for achieving quasi-phase-matching using this scheme are explored and the growth of the harmonic 
intensity as a function of experimental parameters are investigated. 

Please note that this is an arXiv version of the original APS paper. Please cite original paper 
L. Z. Liu, K O'Keeffe, and S. M. Hooker, Phys. Phys. Rev. A 85, 053823 (2012). APS link here: 
[http : //pra . aps . org/ abstract /PRA/v85/i5/e053823 
PACS numbers: 42.55.Vc 42.81.Gs 42.65.Ky 



I. INTRODUCTION 

High harmonic generation (HHG) is a nonlinear pro- 
cess in which odd multiples of a fundamental driving field 
are produced when an intense laser pulse is focused into 
a low density gas. HHG is an attractive source of tem- 
porally and spatially coherent, tuneable light with wave- 
lengths in the XUV and soft X-ray range and has found 
applications in areas such as time resolved measurements 
[iHl], ultrafast holography [ij, or diffractive imaging Q. 
A semi-classical theory of HHG has been developed by 
Corkum [6j and a quantum treatment has been given by 
Lewenstein et al 0- 

Although HHG is an attractive source for a wide range 
of experiments, its adoption in many applications is pre- 
vented by low conversion efficiency, resulting in low sig- 
nal strengths. Typical conversion efficiencies for HHG 
are approximately 10^^ at photon energies around lOOeV 
and 10^^^ for photon energies above IkeV. This low con- 
version efficiency is caused by a phase mismatch between 
the driving laser field and the harmonics generated at 
each point in the generating medium. This results in 
the oscillation of the harmonic intensity with propaga- 
tion distance, with an amplitude of 100%, preventing the 
continuous growth of the harmonic field. The period of 
this oscillation is 2Lc, where L^ = 7r/Afc is the coherence 
length, and Afc is the phase mismatch. Phase mismatch 
arises from neutral gas and plasma dispersion and, if one 
is employed, the waveguide used to guide the driving ra- 
diation ^,[9]]. 

When generating harmonics in a hollow-core waveg- 
uide it is possible to balance the dispersion due to the 
waveguide and free electrons with that of the neutral 
atoms. In this case true phase-matching, Afc = 0, can be 
achieved resulting in quadratic growth of the harmonic 
signal over extended regions. However, this approach is 
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limited to low levels of ionization, and consequently low 
laser intensities, due to the dominance of the free elec- 
tron dispersion at higher ionization levels, placing a limit 
on the maximum harmonic order which can be phase- 
matched using this technique [^, Q . Other phase match- 
ing schemes include difference frequency mixing or angu- 
lar tuning of crossing beams have been proposed |10l . Ill| 

An alternative approach to overcoming the phase mis- 
match is to suppress HHG in out-of-phase regions, a tech- 
nique known as quasi-phase-matching (QPM). Various 
schemes have been developed for quasi-phase matching 
HHG. A series of gas jets, appropriately-spaced, has been 
used to achieve QPM[i2|,[lJ|, but with this technique the 
total number of zones is limited by the Rayleigh range 
of the focused driving laser. Most quasi-phase match- 
ing schemes rely on guiding the driving laser pulse in 
a hollow-core waveguide since this extends the region 
over which harmonics can be generated. Experiments us- 
ing a train of counter-propagating pulses that represses 
or scrambles the harmonic generation in the destructive 
zones have demonstrated quasi-phase-matching over up 
to 5 zones of harmonics with photon energies in the range 
[1Jt13 ■ Other QPM schemes include usin g a corrugated 
waveguide fly] or multi-mode beating [20, [2l[ in hollow- 
core waveguides. 

In this paper we propose a new QPM technique - 
polarization-beating QPM (PBQPM) - which utilizes po- 
larization beating in a birefringent waveguide to modu- 
late the generation of harmonics [22| . The key advantage 
of PBQPM is its simplicity, since it avoids the need for 
additional laser pulses or longitudinally-structured gas 
targets. Instead, the only requirement is a birefringent 
waveguide with a suitable value of the birefringence. 

In Section |TT] the concept of PBQPM is outlined. Sec- 
tion Im] develops the theory of PBQPM, and Section HVl 
explores the optimal parameters of this scheme. In Sec- 
tion|Vl birefrigent wavegudies and advantages of and lim- 
its to PBQPM are discussed. 



II. POLARIZATION BEATING QPM 



III. THE ENVELOPE FUNCTION FOR PBQPM 



It is well known that the single-atom efficiency of HHG 
depends sensitively on the polarization of the driving 
laser field |23l - |26| . which arises from the fact that the 
ionized electron must return to the parent ion in order 
to emit a harmonic photon. 

In PBQPM a birefringent waveguide is used to gener- 
ate beating of the polarization state of a driving linearly- 
polarized driving laser pulse, thereby modulating the har- 
monic generation process. QPM will occur if the period 
of polarization beating is suitably matched to the coher- 
ence length of the harmonics. In a birefringent waveg- 
uide, the incident radiation can be resolved into two 
components polarized along the birefringent axes. For 
linearly-polarized incident light these components are ini- 
tially in phase, but the difference in their phase velocity 
will cause the two components to develop a phase differ- 
ence which increases linearly with propagation distance. 
As a consequence, the resultant polarization state will 
evolve from linearly polarized at an angle O to (for ex- 
ample) the fast axis, through (say) right-handed ellipti- 
cal polarization to linearly polarized at an angle —0 to 
the fast axis, and thence through elliptical polarization 
back to linearly-polarized radiation parallel to the inci- 
dent light. 

The polarization beat length Lb is defined to be the 
distance for the two polarization components to develop 
a phase difference of tt, and hence: 
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where A/3 is the difference in propagation constant for 
the two polarization components, A is the vacuum wave- 
length, and B — -^^ is the dimensionless birefringence 
parameter 

Thus, by matching the beat length to an appropriate 
multiple of ic, harmonic generation can be turned on 
and off along the length of the waveguide in a controlled 
way - enabling quasi-phase matching. Since the harmon- 
ics are generated most efficiently for linear polarization, 
the general condition for lowest-order for PBQPM is that 
the separation, Lh, of adjacent points of efficient har- 
monic generation is equal to an even number of coherence 
lengths, i.e. PBQPM requires, Li, = nLc, where n is an 
even integer As is discussed in more detail below, since 
the coherence lengths for harmonics polarized parallel to 
the X- and y-axes are different, for a given L^ it is only 
possible to quasi-phase-match one of these components. 
As a consequence the output harmonics will be predom- 
inantly linearly polarised along the matched birefringent 
axis. Moreover, as discussed below, Q and the harmonic 
order q determine the width of the generation zone in 
relation to Li, where Q can be optimized. 



If we write the electric field of the qth harmonic as 

Eg{z, t) = S,{z, i)e*[fe(9")^-9-'*] (2) 

then, within the slowly-varying envelope approximation: 



dz 



^j\(;2')g-«[fe(9w)-9fe(w)]z 



(3) 



where ^ is the electric field envelope, k{ijj) = P{uj) is the 
propagation constant for radiation of frequency w, ^ is 
a normalization constant and A(z) is a relative source 
term. 

For each polarization the wave vector mismatch for a 
specific harmonic, q, can be written as: 

Afc = k(quj) - qk{uj) (4) 

and the envelope function can be solved as: 

^{z)^A f\z'e-'^''''A{z'). (6) 

Jo 



A. Effect of polarization state on HHG 

The relative number of harmonic photons generated in 
the gth harmonic by a driving beam of ellipticity s may 
be written as: 
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(7) 



where e is defined as the ratio between the minor axis to 
major axis. 

Within the perturbative regime a — q — I, as verified 
by Budil et al [23] for harmonics q = 11 to 19, and by 
Dietrich et al for harmonics up to g « 31 |23|. Schulze et 
al found that for higher-order harmonics the sensitivity 
of harmonic generation to the ellipticity of the driving ra- 
diation is lower than predicted by Eqn ([7]) with a = q—1 
[28| . although in this non-perturbative regime the effi- 
ciently of harmonic generation still decreases strongly 
with e. Further measurements of the dependence of har- 
monic generation on ellipticity have been provided by 
Sola et al [241 • It is recognized that Eqn ([7]) is an approx- 
imation, but it will serve our purpose of demonstrating 
the operation of PBQPM. 

The offset angle and ellipticity of the harmonics gener- 
ated by elliptically-polarized radiation have been shown 
to depend on the ellipticity and intensity of the driv- 
ing radiation, and on the harmonic order [2^, ^Q, \2ah 
|30| . Propagation effects can also play an important role. 
Since the amplitude with which harmonics are generated 
decreases strongly with increasing ellipticity, we are most 
interested in the ellipticity of the harmonics generated for 



small e. It has been shown that for higher-order harmon- 
ics, and/or high driving intensities, both the ellipticity 
and change in ellipse orientation of the harmoni cs g ener- 
ated by radiation with e ~ are close to zero [2^- We 
will therefore make the simplification that the generated 
harmonics are linearly polarized along the major axis of 
the driving radiation, and that the harmonics polarized 
along the fast and slow axes of the waveguide may be 
treated separately. 



B. Birefringence & evolution of the ellipticity of 
polarization 




In this section, the evolution of the driving field's ellip- 
ticity will be developed. We assume azimuthal symetry 
of the modes. Let x and y be the birefringent axes of 
the waveguide with the x axis being the slower axis such 
that l3x < I3y The driving field can be decomposed into 
the X and y components: 
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(8) 

where E{r) describes the transverse electric field profile 

as a function of the distance r from the propagation axis. 

At a given point z the electric field has an ellipticity 

given by: 



e{z) 



where 
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(10) 



Here (t>max and (j)min are the values oi cj) — tot which give 
respectively the maximum and minimum magnitude of 
the electric field. 

Although the ellipticity is periodic with period Lb, the 
electric field is periodic with period 2L;,. The beating of 
the ellipticity as a function of propagation distance for 
various different angles of incidence is shown in Fig. [TJ 



C. Wave- vector mismatch in a birefringent 
■waveguide 

In a birefringent waveguide the wave vector mis- 
match A/c is polarization-dependent. However, the only 
polarization-dependent term in Eqn (Ois ^kwaveguide, 
the wave vector mis-match arising from the waveguide 
dispersion. Hence we may write. 



^f^x, waveguide 
^f^y, waveguide 
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From this, we may then write, 
Afcx = Afc„ - 



FIG. 1. Ellipticity beating for different incidence angles O. 
From top to bottom: Dotted light gray line shows Q — 45°, 
Solid gray line shows Q — 30° and O = 70°; and dot-dashed 
gray line shows O — 20° and O — 80°. 



where. 
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and the approximation follows if the birefringence is small 
at the frequency of harmonic q. 



D. Constructing the PBQPM envelope function 



From the discussion in Section UlIAl we may treat each 
polarization separately, and hence Eqn (|6]) becomes: 



e.(z) = A/,^dz'e 



ils.kx 



P^x{z') 

Ay{z') 



where 



A.(z0^ y7RiO] sine 
Ay{z') = ^f[e(z')] cos Q 



(15) 



(16) 



Skr 



are the relative source terms. In terms of the faster po- 
larization state (the y component), 

" e. (z) ^ A sin e /„" dz'e-'^^'^y+'i ^^')^' y^ WT] (.j. 
Uz) = AcoseJ^dz'e-'^''y''y/fWY] ^ ' 

Fig. [2] shows an example of /[e(z)] and the real part 
of the envelope function intergrand, assuming a = q — 1 
and Lb = 2Lc^y Notice that integrating across the y 
component would result in a monotonic increase of the 
HHG amplitude whereas for the x component, no net 
increase in harmonic amplitude would result owing to the 
rapid oscillations of the integrand. We conclude that the 
the polarization of the output harmonic will be almost 
perfectly linear, in this case with an electric field parallel 
to the y-axis since it is for that polarization that the 




FIG. 2. Top graph (a): Variation of the harmonic generation 
efficiency, f{z) = f[e{z)], as function z ior q = 27 and Q — 
45°. Bottom graph (b): Real part of the of the x component 
of the integrand (solid gray) and y component (dashed black 
line). Assuming Lb = ^L^y, O = 45° and q — 27 



polarization beating has been matched to the coherence 
length 

Fig. [3] compares PBQPM for these parameters against 
perfect phase matching and perfect QPM; we define the 
latter to correspond to square- wave modulation of the lo- 
cal harmonic generation with a period 2Lc and complete 
suppression of the out of phase zones. 



IV. PROPERTIES OF PBQPM 
A. Coherence length and beat length matching 

For a waveguide with fixed birefringence B the con- 
dition for PBQPM can be realized by tuning the gas 
pressure and/or adjusting the driving laser intensity un- 
til the period of polarization beating and the coherence 
length are appropriately matched. The general condi- 
tion for QPM is that the distance over which harmon- 
ics are generated and that over which generation is sup- 
pressed are both equal to an odd number of coherence 
lengths, i.e. to {21 + l)Lc and {2m + l)Lc respectively, 
where I and m are integers. Thus the QPM period is 
in general 2{l + m + l)Lc = nLc and hence the gen- 
eral condition for PBQPM is Lt, = nLc, where n is 
even. Note that if the distances over which harmonics 
are generated and suppressed are equal, i.e. I = m, then 
PBQPM requires satisfaction of the more restrictive con- 
dition Lb = 2{2l -\- \)Lc. These considerations are illus- 
trated by Fig. Uwhich shows the calculated growth of the 




FIG. 3. Relative HHG intensity for different phase matching 
conditions assuming that Lj, = 2Lc,y, q = 27, Q = 20° for 
PBQPM. Dashed magenta line shows perfect phase match- 
ing, dashed black line for perfect quasi-phase matching, dot- 
dashed cyan line for no phase matching. The solid blue line 
shows the relative amplitude squared for the x component 
whilst the dashed yellow line shows the relative amplitude 
squared for the y component. The thick gray line shows the 
total HHG intensity for PBQPM. 



harmonic intensity for the cases n — 1,2,3 and 4. It may 
be seen that monotonic growth of the harmonics does not 
occur when n is odd, as expected. Notice that PBQPM 
does occur for the case n = 4, which corresponds to the 
lengths of the suppressed and unsuppressed regions of 
harmonic generation being unequal; this is possible for 
PBQPM since the efficiency of harmonic generation is 
very sensitive to the ellipticity, so that it is possible for 
the harmonic generation regions to be shorter than those 
in which generation is suppressed. 



B. Output intensity and coupling angle O 

It can be expected that the coupling angle Q has an 
optimum: if Q is small then the s remains small for all 
z, and the harmonic generation is not sufficiently sup- 
pressed in the out of phase zones; if it is too large then 
harmonics are only generated for a small fraction of the 
in-phase zones. Fig. [S] shows the relative harmonic in- 
tensity as a function of propagation distance for various 
coupling angles 6. Fig. [5^ shows as a function of Q, 
and for various harmonics q, the normalized harmonic 
signal after one beat length. Also shown is the optimum 
coupling angle as a function of harmonic order q. 

The optimum coupling angle, Q, will decrease with 
harmonic order since higher-order harmonics are more 
sensitive to the ellipticity of the driving radiation as dis- 
cussed in Section llllBI . and hence for these harmonics it 
is possible to couple more of the polarization of the inci- 
dent radiation along the polarization of the harmonics to 
be generated by PBQPM. Fig. [SJd confirms this for the 
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FIG. 4. PBQPM for 6 = 34°, g = 27, and Lb = nLc,y with 
n = 1,2,3,4. The top graph shows the relative harmonic 
intensity as a function of propagation distance z for perfect 
QPM (dashed black line). Lb = 4Lc,i, (dot-dashed red line). 
Lb — 3Lc,y (solid green line). Lb = 2Lc,y (dashed blue line), 
and Lb — Lc,y (dotted magenta line). The bottom four graphs 
show the p = y polarization source term Ay as a function 
of propagation distance z for {h)Lb = 4Lc,y,{c)Lb — SLc.y, 
{d)Lb = 2Lc,j;, and {e)Lb = Lc,y 



FIG. 5. PBQPM for q = 27 and Lb = 2Lc,y with 9 = 
30°, 80°, 45°, 20°. The top graph shows the relative harmonic 
intensity as a function of z for perfect QPM (dashed black 
line) and coupling angles of = 30° (dot dashed red line), 
B — 80° (solid green line), Q = 45° (dashed blue line), 
O — 20° (dotted magenta line). The bottom four groups 
show the p — y polarization source term Ay as a function of 
z for (b)e = 30°, (c)e = 80°, (d)e = 45°, and (e) 6 = 20°. 



case of Lh — 2Lc 



V. DISCUSSION 

Having introduced the concept of PBQPM and ex- 
plored its operation in theory, we now discuss how this 



quasi-phase-matching scheme might be realized in prac- 
tice. Several types of birefringent waveguides have been 
developed to date: elliptical or rectangular waveguides 
31 1; waveguides made from birefringent materials |32| : 
waveguides with nonuniform refractive index [3l|; pho- 
tonic crystal fibers (PCFs) |33l - [35| : waveguides con- 
structed with tunable piezoelectric material |36l|: tun- 
able stress-induced birefringent waveguides |3l| . Many 
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FIG. 6. Top graph (a): Relative HHG intensity as a function 
of coupling angle ©, normalized to perfect QPM after one beat 
length where Lb — ^L^y for q = 15 (solid red line), q = 27 
(dashed green line) , g = 39 (dotted blue line). Bottom graph 
(b): Optimal angle as a function of harmonic order q for 



of these birefringent waveguides employ solid cores, and 
are therefore not suitable for PBQPM. Hovifever, we note 
that (non-birefringent) gas-filled PCFs have been used 
to generate high harmonics [37.] . Moreover, commer- 
cially available hollow core birefringent PCFs have been 
demonstrated with birefringence parameters as large as 
B > 10"'' [35|, |3^, corresponding to matching a co- 
herence length of ~ 1mm, and solid-core PCFs with 
B > 10"'^^ have been demonstrated [30|- The birefrin- 
gence of optical crystals can be as large as i? « 0.4 
|40| . Although the birefringence of hollow-core waveg- 
uides made from materials of this type would be only 
a fraction of that of the wall material, it seems likely 
that construction of hollow-core waveguides with suf- 
ficient birefringence to match coherence lengths in the 
range of a few hundred microns will be possible. The de- 
velopment and study of hollow-core birefringent waveg- 
uides of this type will form the basis of future work. 

The principal advantage of PBQPM is the simplicity 
of the experimental set up: a single laser pulse correctly 
coupled into a birefringent waveguide. Other QPM tech- 
niques are likely to be limited by several factors: for 
pulse-train QPM, the number of ultrafast pulses which 
can be generated; the precision with which corrugated 
waveguides or gas jet arrays can be manufactured; or, 
for multi-mode QPM, different damping rates for the two 



waveguide modes. PBQPM does not suffer from these 
difficulties. As for any HHG scheme limits will be set by 
absorption of the harmonics and variation of the ioniza- 
tion level - and hence of the coherence length - within 
the driving pulse; and as for other techniques employing 
a waveguide, damping of the driving radiation may, in 
principle, limit the length over which harmonics will be 
generated. A potential limitation unique to PBQPM is 
dispersion of the two polarization components: the slower 
polarization will lag the faster one until the two polariza- 
tion states separate and polarization beating ceases. In 
principle this limitation could be overcome by reversing 
the birefringence after the two polarization components 
have slipped significantly; the limits imposed by all these 
mechanisms will be explored in future work. 

Finally we note that it is also possible to achieve 
PBQPM in non-birefringent waveguides by exciting 
different-order waveguide modes with nonparallel polar- 
izations. In this arrangement the different phase veloc- 
ities of the modes will lead to polarization beating in 
addition to the beating of the intensity which is utilized 
in muhi-mode QPM ^. Indeed, driving PBQPM with 
different-order waveguide modes increases the available 
parameter space and could allow higher-order harmonics 
to be quasi-phase-matched and greater efficiencies to be 
achieved. 



VI. CONCLUSION 

We have proposed a new method for quasi-phase- 
matching HHG - polarization-beating QPM (PBQPM) - 
in which polarization beating in a birefringent waveguide 
is employed to modulate the efficiency with which har- 
monics are generated. Quasi-phase-matching may then 
be achieved by suitably matching the period of polariza- 
tion beating with the coherence length. 

We have developed a theoretical model of PBQPM and 
used this to demonstrate that PBQPM can indeed en- 
hance the intensity with which harmonics can be gener- 
ated. This model was also used to explore the optimum 
coupling angle for PBQPM and to show that this depends 
on the harmonic order q. 
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